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Abstract 

We present a method to determine tlie stop mixing angle and its CP-violating phase at the 
LHC. As an observable we use ratios of branching ratios for different decay modes of the 
light stop ti to charginos and neutralinos. These observables can have a very strong depen- 
dence on the parameters of the stop sector. We discuss in detail the origin of these effects. 
Using various combinations of the ratios of branching ratios we show that, depending on the 
scenario, one can achieve accuracies in the range of a few percent for determining the light 
stop mass, the mixing angle and the CP phase. 
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1 Introduction 



Supersymmetry (SUSY) [1] is one of the best-motivated extensions of the Standard Model 
(SM). It allows one to stabilize the hierarchy between the electroweak (EW) scale and the 
Planck scale and to naturally explain electroweak symmetry breaking (EWSB) by a radiative 
mechanism. The naturalness of the scale of electroweak symmetry breaking and the Higgs 
mass places a rough upper bound on the superpartner masses of several TeV and the fits to 
the electroweak precision data point to a rather light SUSY spectrum [2]. Therefore a high 
potential for the discovery of SUSY is expected at the LHC [3J. 

Strongly interacting SUSY particles will be produced copiously at the LHC, with cross 
sections up to tens of pb for squarks and gluinos if their mass range is of a few hundreds 
GeV. Cross sections for direct production of scalar top quarks - the supersymmetric partners 
of top quarks - are expected to be smaller due to a different production mechanism, however 
still in a range of a few pb, e.g. for the SPSla' parameter point [3]. The other possible source 
of 3rd generation squarks, depending on the details of particle mass spectrum in the SUSY 
scenario, would be decays of other squarks and gluinos [5]. 

Stops are of a special interest since they play an important role in the mechanism of 
radiative electroweak symmetry breaking. Light stops together with CP-violating phases 
can also provide an attractive mechanism for electroweak baryogenesis by triggering a strong 
first-order electroweak phase transition [6]. Therefore a careful analysis of the stop sector 
can give an insight into the mechanism of EWSB and the origin of the matter-antimatter 
asymmetry. Finally, the stop sector has a large impact on the masses of the Higgs bosons [7], 
and in the presence of the CP- violating phases it triggers mixing between CP-odd and CP- 
even Higgs states [8] . Therefore a precise knowledge of the 3rd generation squark parameters 
would allow us to test the anatomy of the Higgs boson sector and electroweak symmetry 
breaking in the MSSM. 

If stops are within the kinematic reach of the International Linear Collider or CLIC, 
production cross sections can be measured with a high accuracy [HI [121 [13] . Using polarized 
beams, this can provide information on the mixing angle and masses, and a precise deter- 
mination of stop sector parameters can be foreseen [TU [151 HB] • In this paper, however, we 
concentrate on the measurement of the stop sector at the hadron collider. One of the ad- 
vantages of hadron colliders, like the LHC, is the enhancement of cross section due to strong 
interactions [HI [10] • On the other hand, however, further challenges will become important 
due to the harsh experimental environment. 

One of the possible sources of stops at the LHC will be decays of other supersymmetric 
particles, for example gluinos. The analysis of kinematical edges in the invariant mass 
distributions of the cascade decay chains provides an example of one measurement that 
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is often studied at the LHC Taking the SPSla' scenario as an example, a large number 
of stops and sbottoms will appear in the gluino decay chain. Both, however, can give a 
similar experimental signature and consequently one has to do a simultaneous analysis of 
sbottom and stop sectors. This leads to good constraints for the sbottom sector but the 
constraints on the stop mixing angle are much weaker O . Another possible observable 
is the polarization of top quarks in the decay ti Xi t- The information on the stop mixing 
angle can be extracted here from the forward-backward asymmetries in leptonic and hadronic 
top decays |18]. 

In this paper we focus our attention on the decays of the light top squark to charginos 
and neutralinos that are possible in a wide range of scenarios of the Minimal Supersymmetric 
Standard Model (MSSM): 



The stop and sbottom decays have already been analyzed in the literature in some detail, 
including radiative corrections [HI [201 Ell 1221 [23] . In this paper we propose a method to 
measure the properties of the stop sector using simultaneously the decays Eqs. ([T|) and ([2]). 
We analyze three scenarios of the MSSM with different gaugino/higgsino composition of 
charginos and neutralinos. We show that the branching ratios for these decays can be a 
sensitive probe of the mixing angle in the stop sector and also of the CP- violating phase. We 
use a model-independent approach, i.e. without assuming a particular structure for the stop 
mass matrix, and parametrize the stop interactions in terms of the mixing parameters cos 6^ 
and (pi- Since the absolute measurement of branching ratios is expected to be very difficult 
at the LHC we propose to exploit another set of observables - ratios of branching ratios, cf. 
Ref. [SlIlTj. We argue that by looking at direct stop pair production and the following decays 
one can get a good accuracy for the determination of the mass and the mixing parameters 
of stops. We briefly discuss possible experimental issues for these processes. Finally, a fit 
is performed to give a range for the expected parameter determination precision. 

The paper is organized as follows. In Section [2] we give a brief overview of the mixing 
and the couplings of the stop, chargino and neutralino sectors of the MSSM. In Section [3] we 
give analytic expressions for the decay widths of the light stop into charginos and neutralinos 
and analyze their dependence on the stop mixing parameters in chosen scenarios. Section [4] 
explains in detail how to determine the stop mixing parameters using stop decays at the 
LHC for our benchmark models. Finally, we summarize our findings in Section [51 



xtb, 





(2) 
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2 Sparticle mixing and couplings 

2.1 Stop sector of the MSSM 

In the Minimal Supersymmetric Standard Model the stop sector is defined by the mass 
matrix Adj in the basis of gauge eigenstates (f^, tu). The 2x2 mass matrix depends on the 
soft scalar masses Mq and Mu, the supersymmetric higgsino mass parameter /i, and the soft 
SUSY- breaking trihnear coupling At. It is given as 



where 



mL = M^ + mlcos2(5{^-'^sm^9w), (4) 
2 

m^fl = Mir + -m| cos 2/5 sin^ 9^ 



(5) 



rriLR 



= At -fj,* cot P, (6) 



and tan/9 = V2/V1 is the ratio of the vacuum expectation values of the two neutral Higgs 
fields which break the electroweak symmetry. Prom the above parameters only /i and At can 
take complex values 

At=\At\e"^\ l^=\l^\e"^^ (0 < 0^, 0^ < 27r) , (7) 

thus yielding CP violation in the stop sector. 

The hermitian matrix Al| is diagonalized by a unitary matrix TZf that rotates gauge 
eigenstates, tL and tR, into the mass eigenstates ti and 

where we choose the convention mf^ < for the masses of ti and ^2- The matrix TZf rotates 
the gauge eigenstates, II and {r, into the mass eigenstates ii and ^2 as follows 

^1 A ^ f \ ( cos% sin6'£e~"^* \ ( ii 



t2 J~ ^* [tR ) [ - sin Oie'l'i cos Of ) [ tR ) ' 

where 6^ and 0^- are the mixing angle and the CP- violating phase of the stop sector, respec- 
tively. The masses are given by 

^h,2 = I (2771^ + + mlji T \J {ml^ - mlj^y + 4|mM|2mf) , (10) 
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whereas for the mixing angle and the CP phase we have 



a -mt\mLR\ , . 

cosOt = = , (11) 
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ti 



sinOi = (12) 
Jm1\mLR\'^ + (m| - m|^)2 

(t)i = arg(A-/i*cot/?). (13) 

By convention we take < < vr and < 0^ < 27r. It must be noted that 0^ is an 'effective' 
phase and does not directly correspond to the phase of any MSSM parameter. Instead, the 
phase will have contributions from both (pAt and 0^. If ttill < f^RR then cos^ > | and ti 
has a predominantly left gauge character. On the other hand, if > rriRR then cos^ ^? < | 
and ti has a predominantly right gauge character. 

2.2 Chargino mixing 

In the MSSM, the mass matrix of the spin-1/2 partners of the charged gauge and charged 
Higgs bosons, and H^, takes the form 

/ Ms V2mwsmp \ 

Mc=[ ^ . ' (14) 

y v2mvi^cosp /i J 

where M2 is the SU (2) gaugino mass parameter. By reparametrization of the fields, M2 can 
be taken real and positive, while the higgsino mass parameter fi can be complex, see Eq. ((Tj). 
Since the chargino mass matrix Aic is not symmetric, two different unitary matrices are 
needed to diagonalize it 

f rn-± \ 

U*McV^ = with m.± < m^± . 15 

\ m^± / Xi X2 ^ ' 

U and V matrices act on the left- and right-chiral 'ipL,R = {W,H)l r two-component states 

xf = Ujk^^, xi = v,ki^j^, (16) 

giving two mass eigenstates x^, xf. 

2.3 Neutralino mixing 

In the MSSM, the four neutralinos x^ {i = 1,2,3,4) are mixtures of the neutral U{1) and 
SU{2) gauginos, B and W^, and the higgsinos, and H2. The neutralino mass matrix in 
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M 



N 



(17) 



the {B,W^,H^,H^) basis, 

f Ml -mzcpsw mzspsw \ 

M2 mzcpcw -mzspcw 

-mzcpsw mzCfsCw -/i 

\ mzSf^sw -mzspcw -/i / 

is built up by the fundamental SUSY parameters: the U{1) and SU{2) gaugino masses Mi 
and M2, the higgsino mass parameter /i, and tan/5 = t;2/t^i (c/3 = cos/5, svi/ = sin 6^1^ etc.). 
In addition to the fi parameter, a non-trivial CP phase can also be attributed to the Mi 
parameter: 

Ml = \Mi\ e''f", (0 < 01 < 27r) . (18) 

Since the complex matrix A^at is symmetric, one unitary matrix is sufficient to rotate the 
gauge eigenstate basis {B, W^, H^, H2) to the mass eigenstate basis of the Majorana fields Xi* 

diag(m^o, m-^o, m^o, m-^o) = N*AiNN'^ , i^xl ^ "^X2 ^ "^x" ^ "^xli) • 

The masses m^o (i = 1, 2, 3, 4) can be chosen to be real and positive by a suitable definition 
of the unitary matrix A^. 

2.4 Couplings of stops to charginos and neutralinos 

We now give explicit formulae for the couplings relevant for decays Eqs. ([1]) and ([2]) in the 
convention of Ref. [2^]. In terms of two-component (Weyl) gauge eigenstates the coupling 
between stop, top and neutral gauginos/higgsinos is given by 

C = iV2 il (^g2T^W' + ^giB^ II - i'^gii*j,Btn - Yti*^H% - YtilH%^ + h.c. , (20) 

where e = g2Sw = giCw, = is the SU{2) generator and is the Pauli matrix. After 
electroweak symmetry breaking for the mass eigenstates t, ti and we get 



A?xO = X° (QfPL + Q'fPR) t i* + h.c. , (21) 



where 



Q-f = 'Ril(\swN*l + cwN*^ - Y.'RJ.^N*,, (22) 

Q-f = ^ 7^*2iv,l - 7^fliv,•4 , (23) 



with the top Yukawa coupling given by 



= • (24) 

V2 mvp^sv^sin/? 
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We now see that the right squark couples only to the bino and the higgsino components 
of the neutralino. If the fi parameter is much larger than the gaugino mass parameters 
then the light chargino xf and light neutralinos Xi? X2 gauginos with small higgsino 
components. In this case the Yukawa term in Eq. (!23l) is negligible for stop decays into these 
states. On the other hand, as can be seen in Eqs. fl22|) and fl23|) . left squarks couple both 
to the bino and the wino, however with the bino coupling suppressed by a factor 1/3 due to 
hypercharge. Therefore, having a prior knowledge on the composition of neutralinos we can 
infer the structure of the stop sector by comparing strength of the stop coupling to different 
neutralino states. 

Let us turn now to the coupling between chargino, stop and bottom quark. The interac- 
tion Lagrangian in terms of gauge eigenstates reads in Weyl notation 

£ = ig2ilW^bL + Yti*j,H+bL + YdlH^bn + h.c. (25) 

After electroweak symmetry breaking and rotation of fields to their mass eigenstates we get 

^bix = W {QP'^Pl + QP^'Pr) b i* + h.c. , (26) 

where 

Qg'^ = -— 7^*lV-l + l^7^*2V^, (27) 

■' Sw 

Qj'^ = n7^*lf/,2, (28) 

with the bottom Yukawa coupling given by 

n = ^= — - . (29) 

V2 miysvi/ cosp 

The right stop couples only to the higgsino component of the chargino via the Yukawa 
term in Eq. (1281) . whereas the left stop couples both to the higgsino and the wino. When 
the light chargino is mainly wino-like the higgsino couplings are small and only the stop- 
bottom-wino term becomes relevant. Therefore, similarly as for interactions with neutralinos, 
measurement of coupling strength to the light chargino can probe the left-right composition 
of the light stop. 



3 Stop decays to charginos and neutralinos 
3.1 Analytical formulae 

We start this section by giving formulae for the decay widths of the top squarks into charginos 
and neutralinos [IHIES]. Using couplings defined in Sec. I2.4l the tree-level width for the decay 
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Eq. ([T]) can be written as 



r(t. - xjb) 



IGvrm? 



QTr\ +m 



±,R|2 



-4Re 



2 2 2 

rrir — mi, — m-± 



mbm~± 



(30) 



with the kinematic triangle function 

K,{x,y,z) = ^/{x-y- zy - Ayz, 



(31) 



and the couphngs Qfj given by Eqs. Substituting the exphcit matrix elements 

of Eq. ([9]) we can make the following expansion in terms of the stop mixing angle and the 
phase 

o2 



Re 



li I 



:cos2^,- fn'l^.2|' + ^l^.il') + sin2^^,-r,2|^,2| 



(32) 



' ij ^ij 



- 2 sin Oi cos — Re [e'^'^'^VjiV*] , 

- cos^ e-t — n Re [?7;2^/i] + sin 6^ cos Re [e-''^*-[/;2^/2] . (33) 



We see explicitly that the dependence of the phase (pi appears only if there is a significant 
higgsino component (Uj2 or Vj2) in the chargino we are interested in. 
Analogously, for decays to neutralinos we have [151 ES] 



K(m? ,m;, m|o) 
IGvrm? 

'■1 



r^0,L\2 |^0.i?.|2 \ I z 2 z 



-4Re 



rritm^o 

A. i 



(34) 



with K{x,y,z) given by Eq. (13T!) and couplings Q^, by Eqs. (122!) . (!23|) . Similarly we obtain 



cos^ 9 



* * 2s2 r2 



+ r/|iV,4|^ + sin^ -^|iV,-i|^ + F/|iV,-4 



8e2 



2 sin 6i cos 0^ 



Re 



SwCw 



i2 



e'** ( -swN*, + ctyiV;^ ) iV,4 



qp2 



2v^( 
3cvi/ 



Re[e-'**-iV,-iiV;J 



(35) 



Re 



= COS^ 6: 



V2 



SwCw 



;SwN*i + CwN*^ N*^ 



J. ■ 2ZJ 2V2e 

+ sm e-t— 

ocw 



YtRe[N;,N.,, 



2 e' 



+ sin^i-cos^j Re re-'<^*W*|l ^ Re 



(36) 
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An interesting feature of Eqs. (130|) and is the relative importance of the squared 
|Qf^ |^+|(5^|^ terms and the left-right interference Re [QfjQ^*] terms. As they are multiplied 
by mass factors, it is going to be sensitive to the mass splitting between stop and xt Xi* ^ 
pairs. If the given decay mode is close to its kinematic threshold (which will be the case for 
heavier neutralinos) the second term will become dominant, whereas far from the threshold 
the first term will usually be much larger. 

3.2 Discussion of typical mixing scenarios 

In order to analyze the dependence of the stop mixing angle on the decay widths and the 
branching ratios we consider three benchmark points of the MSSM. The first scenario is the 
well known mSUGRA inspired SPSla' parameter point [4J - in the following we will refer to 
it as Scenario A. A feature of mSUGRA scenarios is that the charginos and the neutralinos 
are to a large extent pure gaugino/higgsino states: the lightest neutralino is bino-like, the 
light chargino and the second neutralino are winos, and the heavy chargino and the heavy 
neutralinos are higgsino-like. Scenarios B and C are adopted from Ref. In Scenario 

B the wino mass parameter M2 and the higgsino mass parameter /i are of a similar order, 
giving strong mixing between the wino and the higgsino components of the charginos and 
the neutralinos. This makes the determination of 9i more difficult since both left and right 
couplings of Eqs. (!2T|) and (!26|) come into play for any value of cos 9i. On the other hand 
this gives the possibility to study the dependence on the CP-violating phase 0^, thanks 
to the last terms of Eqs. ( l32l) . ( 133|) . ( l35i) and ( l36l) . Finally, Scenario C features the wino 
mass parameter two times larger than the /z parameter. In this case higgsino-like states 
will be lighter than winos with rather small mixing. In both cases. Scenarios B and C, the 
lightest supersymmetric particle Xi is bino-like. In order to study the possible dependence 
of branching ratios on the CP-violating phase in the last two scenarios we introduce a CP 
phase for the stop trilinear coupling At. For all three scenarios we keep the values of other 
parameters (i.e. slepton and squark sectors) as in the SPSla' scenario. The values of the 
gaugino, higgsino and stop sector parameters are collected in Tab. [Hand the nominal values 
of masses, mixing angles and branching ratios are listed in Tables [2] and O 

We now discuss the behaviour of the decay widths and the branching ratios with respect 
to the stop mixing angle and the CP phase in each of the scenarios. 

Scenario A - mSUGRA 

According to the discussion in Sec. 12.41 for Scenario A we expect that if the ti is mainly 
a left stop (i.e. for cos 6'^ = ±1) then it will dominantly couple to xt ^"^^ X2 (which are 
both winos), whereas the coupling to the bino-like Xi is suppressed. On the other hand. 
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Ml 


M2 


fi 


tan P 


At 


Mq 


Mu 


Scenario A 


103.3 


193.2 


396.0 


10 


-565.1 


471.4 


387.5 


Scenario B 


109.0 


240.0 


230.0 


10 


-610ei^/2 


511.0 


460.0 


Scenario C 


105.0 


400.0 


-190.0 


20 


-610ei^/4 


511.0 


460.0 



Table 1: MSSM parameters of Scenarios A, B and C relevant in the present study. Mass 
parameters and trilinear coupling are given in GeV. 







^2 










X3 




Scenario A 


366.5 


585.5 


183.7 


415.7 


97.7 


183.9 


400.5 


413.9 


Scenario B 


395.5 


609.0 


178.0 


302.9 


101.4 


182.2 


237.8 


303.2 


Scenario C 


396.9 


608.0 


182.9 


419.0 


99.0 


186.6 


199.5 


418.9 



Table 2: Masses (in GeV) of stops, charginos and neutralinos in Scenarios A, B and C 
calculated by SPheno 2.2.3 [26]. 



Parameter 


Scenario A 


Scenario B 


Scenario C 


cos 9^ 


0.56 


0.62 


0.62 


<l>i 





1.53 


0.80 


r(ti) [GeV] 


1.45 


3.25 


6.36 


BR{i, ^ xib) 


73.5% 


60.7% 


63.7 % 


BR{h ^ xt^) 




17.6% 




BR{h ^ x'lt) 


20.1% 


8.8% 


6.4% 


BR{h X^t) 


6.4% 


12.9% 


8.5% 


BR{h ^ xlt) 






21.4% 


a{pp iiil) [pb] 


3.44 


2.27 


2.27 



Table 3: Nominal values of mixing angles in the stop sector and branching ratios for stop 
decays calculated from Eqs. (!30l) and (IMl) for Scenarios A, B and C. In the last row, cross 
sections for stop pair production at the LHC with ^/s = 14 TeV at NLO from Prospino 

2.1 ^o\m\. 

if the ti is predominantly a right stop (i.e. for cos^f = 0) we should observe enhancement 
in the coupling to the LSP and suppression for the decay to the light chargino and second 
neutralino. This general feature can be seen in the upper left panel of Fig. [H where we 
show the dependence of the decay width on the stop mixing angle cos6'j. The minima for 
decays to chargino and X2 somewhat shifted which is the result of the higgsino Yukawa 
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Scenario A 



Scenario A 



1.000 
0.500 



^ 0.100 
^ 0.050 

0.010 
0.005 



0.001 




-1.0 



0.0 



-0.5 



0.0 

COS 67 



Scenario A 



1 



0.5 



1.0 




PQ 




1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



-0.5 0.0 

cos 6-f 

Scenario A 



3 4 
0f 



1.0 




Figure 1: Decay widths (left column) and branching ratios (right column) for ti in Scenario 
A as a function of the stop mixing angle cos^f (upper row) and the stop CP phase (pi (lower 
row). Black, red and blue lines are for xtb^ x\t and final states, respectively. 



contributions from Eqs. ([22]), ([23]), ([27]) and ([28]). On top of that, the decay ti ^ x^t is 
further suppressed by the phase space, since m^o + = 355 GeV is only slightly lower than 
the light stop mass. As one can see, the decay widths change by an order of magnitude or 
more. Therefore they are a sensitive probe of the mixing between left and right stop states. 
The upper right panel of Fig. [1] shows the dependence of the branching ratios on cos 9i that 
exhibit a similar behaviour as the decay widths. 

Although Scenario A does not contain CP phases, we include them here to analyze the 
sensitivity of the decay widths and the branching ratios. The respective plots can be seen 
in the lower row of Fig. [H The most significant change is for the decay to a chargino and a 
bottom quark. This results from the third term of Eq. fl32l) that changes sign when varying 
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Scenario B Scenario B 




-1.0 -0.5 0.0 0.5 1.0 '0 1 2 3 4 5 6 



COS 6-j (pj 

Figure 2: Decay widths for ti in Scenario B as a function of the stop mixing angle cos 6^ (left 
panel) and the stop CP phase 0^ (right panel). Black, red, blue and green lines are for 
Xi X2 ^ ^"^^ xt^ final states, respectively. 

0£ from to vr giving destructive interference. Although the dependence on 0^- is clearly 
visible the constraints on this parameter, as we will see it later, will be rather weak. 

Scenario B — mixed gaugino/higgsino xt ^nd X2 

The situation changes significantly in Scenario B. The second chargino xt is now lighter 
than the stop ti so there is a new decay channel open. Both charginos and the neutralino 
X2 now have a significant higgsino component. The dependence on the stop mixing angle of 
all the decay widths is much flatter now, as can be seen in the left panel of Fig. [21 We note 
a well pronounced enhancement of the decay widths for right stops (around cos6j: = 0). For 
charginos it can be understood by looking at Eq. fl25|) where the coupling of tji is proportional 
to the large top Yukawa coupling, whereas the coupling of is proportional to the smaller 
bottom Yukawa coupling. For the decay to X2^ the enhancement is due to the left-right 
interference term, whereas for x^t it is due to the quadratic term in Eq. flMl) . 

Thanks to the presence of the higgsino component in charginos and neutralinos, we now 
become more sensitive to the phase 0^, see the right panel of Fig. [21 Apart from the decay to 
X2 1, all other decay widths can change by up to an order of magnitude depending on the CP 
phase. The former remains almost unchanged due to the accidental cancellations between 
two terms of Eq. ([M]) . For the decays to charginos Xi" ixf) the suppression (enhancement) 
of the decay width with the phase arises due to change of the sign of coscp^ when (p^ ^ tc, cf. 
Eqs. ( 132|) and (!33|) . The difference for xf and xt is the result of the sign difference between 
VnV* and l^iK*. 
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Scenario C 



Scenario C 




-1.0 -0.5 0.0 0.5 1.0 



COS 




1 2 3 4 5 



Figure 3: Decay widths for ti in Scenario C as a function of the stop mixing angle cos^^ (left 
panel) and the stop CP phase 0^ (right panel). Black, red, blue and green lines are for xfb, 
Xi X2 ^ S'^d X3 i final states, respectively. 

Scenario C — higgsino-like xti X2 ^^^d X3 

The last discussed scenario features the hierarchy Mi < fi < M2. Therefore the light 
chargino and neutralinos X3 ^^e higgsino-like with small mass differences between them. 
The lightest neutralino is bino-like as in the previous scenarios. The dependence of the decay 
widths on the stop mixing angle has been shown in the left panel of Fig. [31 The difference in 
the decay to the chargino for left and right stops is a consequence of the Yt coupling for right 
states in Eq. (!25|) . A similar effect was seen in Scenario B, however it is now more pronounced 
due to the higgsino nature of the light chargino xf. We also observe the interesting exchange 
of the decay widths to heavier neutralinos when the sign of cos 6^ changes. This feature arises 
due to the Y^^ Re [iV*!] term in the second line of Eq. (l36l) that is enhanced both by the 
large top Yukawa coupling and the higgsino nature of the two neutralinos. Since neutralinos 
X2 and X3 have opposite intrinsic CP parities, cf. Ref. [28], the entries in the neutralino 
mixing matrix that correspond to Xs a^^'s purely imaginary. Therefore the contribution has 
an opposite sign in the decay width and hence different behaviour with respect to the sign 
of cos 6^. 

A similar dependence of the decay widths X2 1 and Xs t on the sign of cos (pi can be seen 
in the right panel of Fig. [31 Its origin is the same as in the above discussed case for cos 6^. 
As before the change in the width of the decay to xtb is caused by change in the sign of the 
last term of Eq. ( |32l) with cos 0^, as 0j is varied from to tt. It is interesting to note that 
now the branching ratio for the decay to chargino xt does not show a strong dependence on 
the phase (pi (as opposite to Scenario B). However the dependence of the branching ratios 
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for the decays to neutralinos is still well pronounced. 



4 Analysis at the LHC 
4.1 Ratios of branching ratios 

As one can see in Figs. [1], [2] and [3l the decay widths can change by up to a few orders of 
magnitude depending on the stop mixing angle and the CP phase. In addition, the branching 
ratios are also very sensitive to these parameters. However since the measurement of decay 
widths and branching ratios will be difficult at the LHC we propose to analyze the ratios 
of branching ratios. That means comparing the number of stops decaying to one final state 
with the number of stops decaying to another final state. Having three decay modes possible 
we can define the following ratios of branching ratios for each of the Scenarios A, B and C 

Figure H] shows the above ratios of branching ratios in Scenario A as functions of cos 9i and 
the CP-violating phase 0^ For Scenario B we have three additional combinations due to the 
decay ti xt^ being open 

BR{h ^ xth) _ BR{h^x\t) _ BR{h^xlt) 

""'^ - BR{h - xth) ' ~ BR{h - xth) ' ""'^ ~ BR{h - xth) ' ^''^ 

For Scenario C due to the decay ^1^X3^ being allowed we have 

BRif^^xl^ _ BR{h ^ X\t) ^,,_ BR{h-.xlt) 

""'^ - BR{h ^Xlty ""'^ - BR{h ^Xlty ""'^ - BR{h ^ xlt) ' ^''^ 

Because of the higgsino nature of neutralinos X2 ^"^^ X3 they are very close in mass and it 
might turn out that they are impossible to disentangle at the LHC. Therefore we define two 
additional ratios by combining the decay modes to X2 ^ ^"^^ X3 ^ 

0I6 ^ BR{h ^ xth) 1, ^ BR{h ^ X\t) 

BR{h ^X%t) + BR{h ^xtit)' BR{h ^ xlt) + BR{h ^ xlt) ^ 

In our analysis we focus on direct stop production pp tit\ in order to have better control 
over the number of observed stops and to reduce the background due to bottom squarks. 
In the SPSla' scenario the cross section for this process amounts to 3.44 pb at the next-to- 
leading order [TOl [27], whereas the total SUSY cross section is 60 pb. The cross sections 
for stop pair production in Scenarios B and C are given in Tab. El Due to mass splitting 
between stop states the cross section for pp ^21^2 is much smaller with a value of 0.26 pb. 
Similarly for sbottoms we get a{pp hih\) = 0.6 pb and a{pp — * 62^2) — 0-^ pb. This 
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Scenai-io A Scenario A 




-1.0 -0.5 0.0 0.5 1.0 1 2 3 4 5 6 



COS (p-^ 

Figure 4: Ratios of branching ratios as a function of stop mixing angle cos% and CP phase 
(pl in Scenario A, left and right panel, respectively. Black, red and blue lines are for ratios 
Rft and (see Eq. ([37D), respectively. 



gives a relatively clean environment for the observation of direct light stop pair production. 
Possible final states are as follows: 

VV Ul ^ txl + ixl ttrx\ + itrxl ^ 4£ 4j 26 + E^^, , (41) 

pp ^ U\ ^ txl + txl ^ ttrx\ + tX? ^ 2£ 4j 26 + E^iss , (42) 

^ ^ + ^ i^+r x? + ^ 3£ 2j 26 + E^i,, , (43) 

PP ^ titl ^ tX? + tX? ^ 4j 26 + E^iss , (44) 



pp 
PP 



till -> txl + 6x^ ^ tx? + htval t 2j 26 + , (45) 
iiil ^ 6xr + &X^ ^ ^^"'^a? + &^+^^^X? ^ 2£ 26 + E^iss • (46) 



The production process itself can be tagged using a clean decay mode for one of the 
stops, for instance the decay to followed by a leptonic neutralino decay and hadronic 
top decay. For an integrated luminosity of £ = 100 fb~^ we would have more than 300 000 
stop pair production events. Assuming that on average 10% of charginos and neutralinos 
decay to leptons in our scenarios [26], taking into account the hadronic top branching ratio 
and a selection efficiency of 5%, cf. Ref. [18], one can expect more than 1000 stop events to 
be observed. Therefore in our further analysis we will assume that 1000 events have been 
correctly identified and show that even with this amount of experimental data one can still 
get strong constraints on the stop mixing angle and the mass. 

The other important point we wish to emphasize are the branching ratios for decays 
of the chargino xf and the neutralino X2 into leptons. Although one may expect that 
the related uncertainty will cancel out to some extent in the ratio i?2t (as in our scenarios 
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xf and X2 have similar gaugino/higgsino composition), this is not true for the other ratios 
involving decays to the LSP. Because our focus here is on the stop sector we will assume that 
the leptonic branching ratios of the charginos and neutralinos are known. However as this 
would require better knowledge of the structure of the gaugino/higgsino sectors it is possible 
that the measurements from the LHC would have to be supplemented by the linear collider 
experiment, where charginos and neutralinos can be measured with a high precision. This 
would be an interesting example of LHC/ILC interplay [17], in particular for the scenarios 
where direct stop production is beyond the kinematical reach of the ILC. 

A large number of SUSY and SM backgrounds are expected for stop production at the 
LHC. The most severe Standard Model background, especially for the channels Eqs. (Hlj)- 
fH6|) . will be ti production. As shown in Ref. [18], for the process Eq. fj44|) this background can 
be effectively removed by using appropriate cuts. In any case, the key feature to distinguish 
the signal from SM background will be missing transverse energy, which is much larger for 
stop production due to large energy carried by the LSPs. 

The most important SUSY background process is going to be gluino production with 
subsequent decays to stops or sbottoms. One important difference between the signal and 
these backgrounds is the number of 6-jets. The signal event always results in exactly 2 b- 
jets, whereas SUSY backgrounds will typically have 4 6-jets and this feature can be used to 
suppress them. 

Finally, we note that the signal process with leptonic top decay, e.g. Xi^ ^ &^ + -E'miss, can 
give the same final state as the decay mode with charginos, i.e. Xib ^ M + -Emiss- However, 
we note that this complication does not affect the result of the fit since it does not introduce 
any new unknown parameters. The fitted observables would be a linear combination of the 
original ones, Eq. (|371) . and the fit would rely on the same set of information. Hence, one 
can combine the above channels and actually enhance the signal. 

An important note is that it will not be sufficient to simply remove as much background 
as possible using the relevant cuts. We will also need to understand with a high degree 
of accuracy how each individual signal channel will be affected by the backgrounds. Un- 
derstanding the background well is required, as for each channel we study, the number of 
background events contaminating the sample will be different. Therefore the pollution due 
to backgrounds will affect our ratios of branching ratios. The reconstruction efficiency, cuts 
and triggers will also have a different effect on each channel and will have to be well under- 
stood for our measurements to be accurate. We leave detailed analysis of these effects for 
our different final states and the additional uncertainties they may induce for future work. 
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4.2 Determination of stop mass and mixing angle 

In order to show the possible advantages of using ratios of branching ratios for the analysis 
of the stop sector we first define the normalized ratios 

^5 = i?} - i?f """'""^ with i?f°'"'°"i = i?}(cos^?°'"'"">), (47) 



where is the actual mixing angle in the given scenario and i, j run over all possible 



channels, i.e. lb, It, 2t etc., cf. Eqs. (EZD-dlOl). According to this definition i?i(cos0?°™''^') = 
0. Furthermore, we assume that we have n = 1000 of well identified events of stop tit^ pair 
production. We now take the expected number of events in each decay mode Ui = n x BRi. 
Note that 1 = ^ only if decays to charginos and neutralinos are the only possible decay 
channels. However, for our method it is not necessary to measure all possible decay modes. 
The statistical error for is Astat'^i = '^^^ resulting error for ratios of branching ratios 
is given by 

Am = J(^]\(^y with m = !^. (48) 



"3 

Before analyzing the expected accuracy of determination of stop sector parameters let us 
study the possible influence of the gaugino/higgsino sector parameters, taking as an example 
Scenario A. The precise knowledge of the LSP mass and the mixing angles of the charginos 
and the neutralinos may only be accessible after the results from a linear e^e~ collider are 
available. In Fig. Owe show the dependence of the normalized ratios Eq. (1371) on the gaugino 
mass parameter M2 and the mass of the LSP, mLSP = f^x^- In the first case we keep the mass 
differences m^o — m^o and m-± — m^o fixed as these are expected to be measured with high 
precision at the LHC. As can be seen, the value of R^^ is very stable in both cases, whilst R^^^ 
and i?2t exhibit an increase for larger values of M2 and mLSp. This is because both i?2t 
i?2t include a branching ratio for the decay to x^*^ that is close to its kinematic threshold. 
Therefore, for an increasing rriLSP or M2 (note that M2 ~ m^o in Scenario A) we approach 
the point where this decay becomes impossible. High sensitivity of the decay width near the 
threshold means that to use such a decay mode to determine the mixing angle, one would 
have to know the masses extremely precisely. In this case the ratio of branching ratios is no 
longer a good observable. Moreover the branching ratio for such a decay usually becomes 
very small. 

In order to analyze the possible accuracy in extracting the mixing parameters of the stop 
sector we start with the example of Scenario A. In Fig. [6] we show the behaviour of the 
normalized ratios of branching ratios, Eq. (H7|l . near the nominal value of the mixing angle 
cos^^?"™™'^^ Using only one of three possible ratios, the smallest error and hence the best 
estimate we get is using the ratio which depends on the dominant decay modes x\ t and 
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Figure 5: Dependence of the normalized ratios of branching ratios -Rj, Eq. (147|) . in stop 
decays on the mass of the hghtest neutrahno (left panel) and the gaugino mass parameter 
M2 (right panel) in Scenario A. The values of the other parameters and the mass difference 
between the LSP and X2? X\ fo^' the left plot are fixed to their nominal values. Red lines 
are for black lines for i?2t! ^^^d blue lines for i?2t) see Eq. (ITTI) . The upper and lower 
lines show the band due to \-a uncertainty with respect to the central value as defined in 
Eq. (gHD. 
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Figure 6: Normalized ratios i?*, Eq. (H7|) . near the nominal value of cos^^ = 0.56 for Sce- 
nario A. Red lines are for black lines for i?2t) blue lines for i?2t! see Eq. f l37l) . The 
upper and lower lines show the l-cx uncertainty as defined in the text. The region allowed 
by all three ratios is shaded in grey. 



Xa^- For the ratio -Rgt ^^e impact of the error is slightly larger due to the limited statistics. 
On the other hand the ratio i?2t gives the weakest constraints because both xl and X2 
winos, hence their couplings to ti follow a similar pattern. We assume here that the values 
of the other SUSY parameters, including the t\ mass, are known. Using only the information 
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Figure 7: A fit of ratios -Rj^, -Rgt ^"^^ -^l- f|37|) . to the stop mass and the mixing angle 
cos^£ (left panel), and the stop mixing angle and the CP phase (p^ (right panel) in Scenario 
A for n = 1000 events. Bold, normal and dashed lines are for 1-, 2- and 3-cr contours, 
respectively. 



from the ratio we get the estimate 

0.53 < cos 0.59 (49) 

at the 1-cT level, cf. Fig. El 

Having at hand three possible decay modes we can constrain not only the mixing angle 
cos % but also the mass of the light stop quark and the CP- violating phase 0^. This can be 
done using a fit defined as follows 

{hj} \ ^ / 

where the error is defined by Eq. (148|) and the sum runs over the respective ratios for each of 
the scenarios, e.g. {i,j} = {{lb, It}, {lb,2t}, {lt,2t}} in Scenario A, cf. Eqs ([37]) -(iO]). The 
results of fitting the stop mass m^^ and the mixing angle cos 9]; in Scenario A are shown in 
the left panel of Fig. [71 We find two minima of that fit the input data well. In order to 
resolve the two-fold ambiguity, additional observables will be needed. Assuming that we can 
pin down the correct solution we get the following 1-a estimate of the two parameters 

mi^=366tl, cos^t = 0.56 ±0.04, = 0.98 ± 0.05 (51) 

for 1000 events. The better lower bound for the measured mass is a consequence of the 
earlier discussed small difference between m^^ and m^o + rrit. In the right panel of Fig. [7| 
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Figure 8: A fit of ratios R^, R^ and Eq. ([37D and Eq. ([38]), to the 

stop mass and the mixing angle cos^^ (left panel), and the stop mixing angle and the CP 
phase (right panel) in Scenario B for n = 1000 events. Bold, normal and dashed lines are 
for 1-, 2- and 3-0" contours, respectively. 



we show the results of the fit to the mixing angle and the phase 0^. As expected, the 
sensitivity to the CP phase is poor and taking into account the possible ambiguity in the 
mixing angle cos 6i, the full range of phases remains allowed. 

The situation changes for Scenarios B and C. We now have 6 possible ratios in each case, 
for Scenario B: R\l R^ R^, R^, R^, and for Scenario C: R^ R^ R^, R^ RH, 
i?3*. The results of the fit for n = 1000 events have been shown in Figs. [H]and[9l We again 
consider two cases: fitting of the mass together with the mixing angle cosd^ and fitting 
of the mixing angle together with the CP-violating phase. In each case we assume that the 
value of the third parameter is known. Charginos and neutralinos now have a significant 
higgsino component and, as we saw in Figs. [2] and [31 the dependence on the mixing angle 
is much weaker. Therefore the constraints for the mixing angle and the mass that we get 
are not as good as in the case of Scenario A. It is interesting to note that in general the 
results of the fit are better in Scenarios A and C (gaugino and higgsino, respectively) than 
in Scenario B (mixed case). Consequently we conclude that the scenario with strong mixing 
between gauginos and higgsinos would be the most difficult to resolve. 

Analyzing both the mixing angle and the phase, we obtain four allowed regions. Never- 
theless only small regions are allowed for the CP phase as our observables are more sensitive 
to it than in Scenario A. Branching ratios are CP-even observables, therefore they cannot 
resolve ambiguities for the CP phase. This shows that for precise measurements in the stop 
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Figure 9: A fit of ratios Rl^, R^ and Eq. ([37D and Eq. ([39D, to the 

stop mass and the mixing angle cos^^ (left panel) and the stop mixing angle and the CP 
phase (right panel) in Scenario C for n = 1000 events. Bold, normal and dashed lines are 
for 1-, 2- and 3-a contours, respectively. 



sector one has to use other CP-sensitive observables, like triple products of momenta [25|[29]. 
Only such a combined analysis of CP-even and CP-odd observables can give an unambiguous 
determination of the stop sector parameters. 

Finally we discuss the results in Scenario C when combining decay modes X2 ^ X3 
as the two close-in-mass higgsino-like neutralinos may be difficult to disentangle at the LHC. 
The fit is now performed to 5 ratios of branching ratios: i?2t, R^l^ Rllt Rllty Eqs. ( 1371) 
and (HOl) . In Fig. [10] we show the fit to the stop mass and the mixing angle (left panel), and 
to the mixing angle and the CP phase (right panel). It turns out that we lose sensitivity to 
the elements of the stop mixing matrix. In such a case additional input, for example from 
linear collider, would be needed in order to resolve properties of the stop sector. 

5 Conclusions 

In this paper we have analyzed the stop sector of the Minimal Supersymmetric Standard 
Model. In particular, we have looked at the couplings and the decays of the supersymmetric 
top partners to the charginos and the neutralinos. As stops play an important role in the 
MSSM, it is crucial to measure their couplings and masses at future colliders in order to 
understand the underlying model. Therefore we have proposed a promising method for the 
determination of the stop sector parameters at the LHC. 
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Figure 10: A fit of ratios E)^, RH R^^ and i?^*^, Eq. ([37D and Eq. (gO]), to the stop 
mass and the mixing angle cos6'f (left panel) and the stop mixing angle and the CP phase 
04 (right panel) in Scenario C for n = 1000 events. Bold, normal and dashed lines are for 1-, 
2- and 3-cr contours, respectively. 



A careful analysis of the couplings of scalar tops to electroweak gauginos and higgsinos 
shows a strong dependence on the mixing angle and the CP-violating phase of the stop sector. 
This effect arises due to the structure of the electroweak gauge couplings and the Yukawa 
couplings of left and right stop states. We have analyzed three benchmark scenarios with 
different structures for the gaugino and higgsino sectors, where the light charginos and neu- 
tralinos had gaugino-like, higgsino-like or mixed composition. Analysis of the decay widths 
and the branching ratios has shown a strong relation between the stop mixing parameters 
and the decay pattern in each of the scenarios. 

Next, we have discussed a possible approach to determine the light stop mass, the mixing 
angle and the CP-violating phase at the Large Hadron Collider. Since stops will be produced 
in large numbers at this machine one can hope to learn the stop properties from their decay 
pattern. As the branching ratios are going to be difficult to be measured at the LHC we 
propose to analyze the ratios of branching ratios for different decay modes. These observables 
inherit a strong dependence on the mixing angle from stop decay widths and therefore can 
be a sensitive probe of the stop sector. Since they rely only on the relative numbers of stops 
decaying via various channels, many experimental uncertainties will cancel. In particular, 
one does not need to control all of the possible decay modes. In fact, as we have shown for 
the SPSla' parameter point, using only two decay modes can give good constraints on the 
stop mixing angle. Finally we have performed fits to show that the ratios of branching 
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ratios can give strong bounds on the parameters of the stop sector: the mass of ti, the stop 
mixing angle cos 9^ and the CP- violating phase 0^. The expected accuracy depends upon the 
scenario studied but looks the most promising for mSUGRA models. 

Application of this method will require the study of many possible final states. Therefore 
a good control of detector effects, like fake rates for leptons and 6-jets, and SM as well 
as SUSY backgrounds will be needed. It is clear that more detailed experimental studies 
are needed to assess the validity of the method and its possible accuracy. However, taking 
into account the importance of the stop sector for our understanding of the supersymmetric 
model, we think that such a study deserves further attention. 
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